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Abstract 

We analyse quasiprobability distributions related to the discrete Heisenberg-Weyl 
f-*) . group. In particular, we discuss the relation between the Discrete Wigner and Q- 

' functons. 
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§ 1 Introduction 

Due to the recent development of hardware facilities, the description of the state and evolution 
^ ! of the wave processes has also evolved. In addition to the configuration space description 

^ I and Fourier analysis, it has now become possible to use the Quasiprobability distributions 

CT' that describe the state of the wave process in phase space of the corresponding mechanical 

system. By wave process we mean any acoustical, optical or quantum mechanical process, 
^ , where, correspondingly, time and frequency, the coordinate of the point where the ray of light 

^ ! intersects the screen and the ray direction, the coordinate and momentum of the quantum 

mechanical particle play the role the phase space coordinates (see, e.g. [1-6]). Until quite 
recently, the practical use of the Quasiprobability distributions was strongly limited by the 
speed and memory restrictions of the computers, but now the situation is changing. The 
possibihty to use the Quasiprobability distributions is highly attractive for it allows us to 
visualize the process in a intuitively clear way: to analyze the signal {e.g. to discriminate 
between its different components, possibly originated from different sources) and to process 
the signal (say, removing noise) at the level which is much higher than the one that can be 
achieved when working with the signal itself or its Fourier transform. 

The theory of Quasiprobability distributions in its most complete form was developed 
in the framework of quantum mechanics {of. [1-4,6]). In its original version, it refers to 
continuous variables and is intrinsically related to the Heisenberg-Weyl group of translations 
of phase space. However, any numerical work involves the discretization, which by itself is 
not a trivial procedure, as is clear from the example of the approximation of the Fourier 
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integral by the Finite Fourier transform (see, e.g. [7]). Therefore, it is desirable to develop 
a systematic theory of the Quasiprobability distributions in discrete phase spaces. The basis 
of such a theory was given in the work of Woottcrs, [8] (see also [9], [10], [11], [12], [13]). 
The Hilbert space of system states was chosen to be a space of periodic functions with a 
finite number of Fourier harmonics. This is obviously a finite-dimensional space. This type 
of quasiprobability distributions is naturally related to the discrete Heisenberg-Weyl group 
of translations of discrete phase space. (It should not be confused with the quasiprobability 
distributions for spin systems, which arc related to the SU(2) group [3, 14]; in this case, the 
finite-dimensional Hilbert space includes spherical harmonics on the sphere.) 

Our goal in this work is to establish the correspondence between the two most important 
discrete Quasiprobability distributions: the Wigner-Wootters ( "W- function" ) distribution 
and the time dependent spectrum ( "Q- function" ) , as it exists for the continuous case, where 
the Q-function can be produced by smoothing the W-function with an appropriate Gaussian. 
A natural periodic counterpart of the Gaussian is the Jacobi Theta-function and therefore, 
they appear in our construction. 

The paper is organized as follows. After the presentation of the discrete coherent states 
(which naturally involves the discrete Jacobi Theta-functions) and the Q-function, we define 
the discrete Wigner function (which is made in terms of the discrete Heisenberg-Wyel group 
and is complementary to the original work [8]). Finally, we establish the relation between 
them. 

We would like to note that we do have in mind possible applications to signal processing. 
However, at this stage of theoretical development, we found it reasonable to confine ourselves 
within the language of Quantum mechanics. 



2 The model 

We consider periodic functions f{x + L) = f{x). Therefore, the Fourier series coefficients 
contain the frequencies, ujk = fc27r/L, where k is an integer. 

Let us consider functions for which Fourier series contain only M coefficients different 
from zero: k = 0, 1, . . . , M — 1. Then all the information is stored in the values of the 
function in M points on the circle: 

Lm 

2^m = -^, m = 0, 1,...,M-1; 

M-l 

/(m) ^ f{xj = E /,e^^-'=-/^. 

fe=0 

The values of the function f{x) at arbitrary points can be recovered by the samphng theorem. 
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We will use the quantum mechanics notation, introducing the basis |mo) = 



M-1 



I/) = Yl HHf), Hf) = f{xj = f{m). 

m=0 



One can use the matrix representation with 



I/) 



/o 
/i 

/m-1 



The Discrete Fourier Transform of the vector /(m) gives us the wave function in momen- 
tum representation: 



M-1 



1 



- ' M 



m=0 



Using the discrete orthogonality relation (1), 



M-1 



(1) 



m=0 



one finds the DFT coefficients, 



Thus, the periodicity also appears in the momentum. 

It is natural to introduce one-step translations in coordinate. 



Txf{Xm)=f{Xm-l): = CXp ( -i — p^, ) , ^ -id^, = 1), 



and in momentum. 



27r 



Tp^exp i — X , xf{x)=xf{x), Tpf{xm)=r"'f{xm), r = e 



i2%/M 



It is clear, that 

= v^f'(l, 1) = rf,%, T(l, 1) = exp 



L 



Note, that in the course of the discretization of the initial continuous model, L completely 
disappears from all of the formulas. 
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3 Vacuum state 

The periodic analogue to Gaussian functions is the Jacobi 6-function defined as follows (see 
[15] and Appendix A), 

n=-oo V ^ k=-oo I ^ ) 

(Here and below we often use the Poisson transformation.) Generally speaking, Re// > but 
we consider only real values of /x. Recall, that 63(2; + n, ji) = 63(2;, /i), and hence x = zL/n. 

Considering the 0-function on a discrete set of points, = mL/M and = rmr/M, 
we come to a discrete periodic Gaussian, 

e(m) = He,) = e3(^ £ e-^^--"/^-'^"^ 

^ ' n=—oo 

This function will play a role of the vacuum squeezed state (with squeezing parameter /i). 
The choice 

IJ, = tt/M 

corresponds to the vacuum Coherent State. Thus, <^ 1, /xM^ 3> 1 will be of interest. Using 
the discrete orthogonality relation (1) one can check the normalization, 

M-l 

Af = (e,|e,)= |e(0^|/.)|2 = M$:$:e-'^"'-^("-'=^)' 

m=0 n k 

[^3(2/x)^3(2/iM2) + ^2(2/x)^2(2/xM2)] , M = odd 
^3(2/i)^3(/xMV2), M = even 

I2nj [^3(2/i)^^3(/iV2) + ^2(2/.)^4(/^72)]/2, M = odd . . 

\ ^3(2/x)^3(2/x'), M = even ^""^ 



= M 



Here 



for unsqueezed vacuum CS, /i = /i' = n/M. We will assume, that /i — iiqTt/M, //q ~ 1, and 
thus, /i, /i' -C 1 simultaneously with iiM"^ , jj! M"^ ^ 1. 
In Eq. (3) we use the notation, 

eM^QMli), A; = 1,2, 3, 4; 
(see Appendix A). For instance. 



In passing from the second to the third hne of Eq. (3) we again used the Poisson transfor- 
mation, 9s{iiM'^/2) = ^2/^777^3 (2//'), and O^i^iiM'^) = ^ n' /2t: 9i{ii' /2). If M > 1, the 
normahzation constant has a simple asymptotic form. Indeed, in this case only the terms 
A; = are important in the Poisson-transformed expressions for ^2,3,4, and neglecting the 
terms O (e-'^^^'/^j, O (e-'^'^'l'^), we have. 



(e>Mie>M) 



2ix 



(4) 



The wave function of the vacuum state in the momentum representation is given by the 
Discrete Fourier transform: 



fe=— 00 



-M(p-feM)2 



TT 



E 



il'K'pnjM—ii'n'^ 



e 



7rp 



4 Discrete Heisenberg-Weyl group 

Let us return to the one-step translations in coordinate and momentum, which are defined 
in the framework of the initial continuous model as — e^^P^^^/^, Tp — ^ix2-k/l ^ 

f^\m) = |m 1), f^f{m) = {m\f^\f) = /(m - 1), f^\p) = r"^|p), 

T,|m)=r'"|m), fpb) = b+l), f,/(p) = /(p - 1). 

Here, once again, r = e*^'^/^. These operators generate a discrete group. It is clear that there 
is a periodicity in momentum, as well as in coordinate, — — 1. These translations 
do not commute, 

rpnrpm nmrpmrpn 

p X X p ' 

Therefore, an arbitrary group element must include the multiplication by a phase factor, 

gis,m,n)=r'f^f;. (5) 

The naive way to introduce the displacement operator in phase space is (compare [10], 
[11]), 

t{m, n) = r'^rn/2frnfn ^ ^-nm/2fnfrn ^ J- _ ^^^^ _ 

However, such an operator does not belong to the group (5). Fortunately, there exists a way 
to improve the situation. Let us now consider the case of odd M. Then, in the set of numbers 
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{n {modM)} = {1, 2, . . . , M — 1} there exists a unique solution to the equation 2n — m, 
which we denote as 



m 
~2 



f m = 2k, 

m = 2k + l. 



(More generally, in the case of prime M there always exists a unique solution to the equation 
sn = m, and the numbers {n (modM)} then form a field.) It is clear, that ^ = y n = 



m and 
We now introduce t 

T{m, n 



+ 

le displacement operator as follows, 

^[nm/2] jimjin _ ^-[nm/2] jm jmi _ ( — ^ 



(6) 



T is unitary, T\m,n) = T ^{m,n) = T{—m, —n), and periodic, T{m + M,n) = T{ni,n 
M) — T{m, n). It is easy to check the multiphcation formula. 



T(m2, n2)T(mi, m) = rl^'"!"^-™)/^]^^^^^ ^ ^ 
One can introduce the adjoint action. 



(7) 



f{q,p)f{m,n)f\q,p) = rP"'-i''f{m,n). 
The matrix elements in the coordinate basis, 

{k\f{m,n)\l) = r"'+[""^/2]4,;+m(WM). 
These matrix elements are orthogonal. 



M-l 



Finally, 



^ {a\T{m,n)\b){d\T{m,n)\c) = M5a4hc. 

m,n=0 



Tr f(m,n) = M6mflSnfl- 



(8) 



(9) 



5 Discrete CS and Q-dunction 

One can generate the complete set of (squeezed) CS by the action of the HW group (6) to 
the vacuum state (2), 



{m\mo,no 



^'j — gi27rno(m-mo/2)/M 



m — mo 1 



^i2nno {m—mo /2)/M 



J2 exp 



t — [mo — m)n — jin 



(10) 
(11) 
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and in the momentum representation, 

oo 

fc=— oo 

The normahzation constant (mo, no, /i|mo, no, /x) — M is the same as in (3). 
The completeness relation holds (for any /x) : 

1 

J? \mo,no,n){mo,no,iJ,\ ^i. 

m(),n()=0 



(12) 



Scalar products of Coherent States. It is enough to consider the vacuum matrix element 
of the displacement operator. 



(e^|mo,no,^) = ( e^|r(mo,no)|e^). 



If M is even, it is equal to 



where. 



{%\m„n,,,) = Mf-^e^ (^1 2^) 6, (^| 2,') , 



a — 3, no — even; j P — 3, mo — even; 
a — 2, no — odd; 1/3 = 2, mo — odd. 



If M is odd, no is even, 
(©^|mo,no,//) = 63 
and if M is odd, no is odd, 
(6/, I mo, no, /x) = 
Asymptotic for large M is 



/Trmo 



-1)™° G3 



Trno 



2M 



2// 64 



Trno 



/Trmo 



V M 



2/x 64 



vrno 



2M 



2" 



2/x Ga 



2M 



vrno 



2M 





mo,no,/i) 







exp 



//'mo + Atno 



Q-function. It is natural to introduce the Q-function as the diagonal matrix element 
of the density matrix between the coherent states. 



Q{m,p) = {m,p,ii\p\m,p,iJ,), 



(13) 



or, for a pure state, p — |^)(^|, 



Q{m,p) = \{'^\m,p,p)f . 
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6 Wigner function 

By analogy with the continuous case, let us introduce the Wigner operator as a two-dimensional 
discrete Fourier transform of the displacement operator, 

■1 M-l 

^(5'P) = T7 E q,P = 0,l,...,M-l. (14) 

-'^ m,n=0 

Therefore, it is an Hermitian operator valued function on the phase space, W{q,p) = W''{q,p). 
One can notice that, 

W{q,p) ^f{q,p)W{0,0)f''{q,p). 
One immediately calculates its matrix elements, 

( i27i 1 
{k\W{q,p)\l) = expj — p(A; - l) \ 52q=k+i(modM)- 

These are precisely the matrix elements of the " Phase Point Operators" considered by Woot- 
ters [8], with the property 

TV {Wiq,p)Wipi, q^)) = M5,,,,5p,p,. (15) 

This property was used to define these operators in [8]. 

The Wigner function (quasiprobability distribution) on the discrete phase space for the 
state I/) is the mean value of the Wigner operator in this state, 

Wf{p,q) = {f\W{q,p)\f). 

To an arbitrary Hermitian operator A, corresponds the Wigner function 

Wf{p,q)^Tr{W{q,p)A). 

If A ^ p is a density matrix and Tr p = 1, then 

J2W{q,p)^l. 

p,q 

For a pure state, Tr = 1, 

Y,W'{q,p) = l. 
p,g 

The orthogonality of the Wigner operator at different phase points (15), has a simple physical 
sense: the Wigner function of the Wigner operator itself is a (5-function, 
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Covariance under HW group. Using the adjoint action of the discrete HW group, (8), 
one can show that 

W{q,p)^f{p,q)W{0,0)f\p,q). 

The Wigner operator at the origin, Vr(0, 0) has a simple form, which we will explain, con- 
sidering examples of odd and even M. For M — 5, 



W{0,0) 



1 

1 

1 

10 

10 



thus the spectrum of Wigner operator for odd M consists of M/2 eigenvalues —1 and M/2 + 1 
eigenvalues +1. For M = 4, 

10 
1 

10 
10 



W{0,0) 



the spectrum for even M consists of (M — l)/2 eigenvalues —1 and (M — l)/2-|- 1 eigenvalues 
+ 1. Therefore, 



M-l 



TrW{q,p) = {m\W{q,p)\m) = 1, 

m=0 



1, 

2, 



M = odd, 
M — even. 



As in the continuous case, W{0, 0) is the inversion operator with respect to the origin. In the 
case of odd M, there is only one stable point (m = 0) under this inversion, while for even M 
there are two stable points m = 0, and m = M/2. In the rest of the paper wc will restrict 
ourselves with the case when M is odd. (In the opposite case the information is lost and one 
cannot reconstruct the state from its Wigner function.) 

Overlap relation. Let us consider the following operator acting in the tensor product of 
two spaces. 



M 



J2W{q,p) X W{q,p). 



p,<i 



Using the orthogonality of the displacement operator matrix elements (9), one can prove that 
a is an exchange operator. 

From here the overlap relation follows directly. Let Wa{p, q) and Wb{p, q) be the Wigner 
symbol of the operators A and B. Then 



— J2WA{p,q)WB{p,q) ^Tr AB. 
p,<i 



(16) 
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Reconstruction of the initial state. If we know the Wigner function Wp{p,a), then the 
density matrix of this state can be reconstructed as 



7 Relation between W and Q-functions 

We assume now that M is odd, and 11 = 11' = tt/M. Let us start with the kernel 

{k\Q(q,p)\l) = {k\q,p){q,p\l) = {k\f(q,p)\e,){e,\f{q,py\l), 
and find its Fourier transform, 



M-l 



^ r""'(A;|Q(g,p)|0 = M^{k\T{m,n)\l)F{m,n), 

q,p=0 



where the function F{m,n) is defined as 




F(m, n) = 
In other words. 



e, 



nm 



TT 



2 



(17) 



M-l 



Q{q,p)= T {m,n)F{m,ny = T{m,n)f{m,ny 



m,m=0 



m,m=— 00 



To find the function /(m, n), let us note that for any periodic discrete function 0(m+M) 
(m) the following formulas hold: 



/ M-l 



TT 



Y (-l)>He- 



m=0 

7 M-l 



■fj,n 



Y{-irHm)e2{nm/M\fj,'), 



n=— 00 
00 



Y 0(^)e 



-M(n-M/2)2 _ / 



m=0 

7m-i 



^ 5^ 0(m)e4(7rm/^l/^'), 



TT 



m=0 



and moreover: 



Y mn])e 



.1 M-l 



—iin 



Y '/•Me>3 



m=0 



TT 
M 



m 
2", 
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/ M-l 



TT 



m=0 



TT 

M 



2". 



Therefore, for any periodic function of two variables, 4>{m + M, n) = 0(m, n + M) = (f){m, n), 
one can prove. 



M-l 



E (p{m,n)F{m,n) = E (p{m,n)f{m,n), 

m,n=0 m,n=— 00 



with 



TT 



fim,n) = (-l)'""exp [-^m'-^n 



(19) 



From here it follows that F{m,.n) — F{n,m) (provided that — /j,', i.e. the coordinate and 
momentum have the same rights. 



8 Conclusions 

The principal results of this paper are the relations (18), (17), (19) between the discrete Q- 
function (time-dependent spectrum) and the discrete Wigner function. As in the continuous 
case, the Q-function can be produced by smoothing the W-function by integrating it with 
the discrete Gaussian function. 

S.C. is grateful to the Department of Computer Science, University of Houston, Texas for 
hospitality. 



Appendix A: Theta functions 

Jacobi Theta functions are the integral functions defined as follows [15]: 
63(2:1//) = 2^ exp|-z2n2; - /in I = J- 2^ exp ' 



n=-oo V ^ fc=-oo I 

II £ f_(.-x(^+l/2)) 

n=-oo V ^ A;=-oo I 

00 ( / -Is /i\2"^ rz" 00 



04{z\ijl) = E (— l)'*exp |— i2n2; — //n^l = w— E ^^P 



2 



e.(.,.) ^ J^exp|-.2(. + l).-,(.+l) £j-l)^exp{- 
ei(z|/x) = f; (-l)-exp|-i2(n + ^)^-/x(n + ^ 



2" 

2. 



E (-l)'exp 
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Here ifj, is quasiperiod, Re /i > 0. Note, that in [15] Theta functions are written as Qk{z, q), 
where r = e~^. ©1,2,4 are shifts of ©3 : 

©4(^|/i) = ©3(^ + vr/2|/i), 02(^|/i) = -^e^"-'^/^©3(^ + ^///2|/x), 
e^{z\ij) = ©2(^-7r/2|/x) = -ie^^-''/^©3(^ + 7r/2 + i/x/2|/^). 

Note, that 

01,2(^ + Vr|/i) = -0i,2(^|/i), 03,4(^ + Vr|/i) = 03,4(^|Ai), 

and 

©1,4(2; + inlii) = -gOi^^izln), ©2,3(^ + ^/^|/^) = gQ2,3{z\n), 
where g — e^'^^+z*. Finally, 

Ql{-Z\l^) = -0i(^|/x), 02,3,4(-^|/^) = 02,3,4(^1/^). 
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